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ABSTRACT 

From the radial velocities of the N iv A4058 and He II A4686 emission lines, and the 
N v A4604-20 absorption lines, determined in digital spectra, we report the discovery 
that the X-ray bright emission line star Wack 2134 (= WR 21a) is a spectroscopic 
binary system with an orbital period of 31.673±0.002 days. With this period, the N IV 
and He n emission and N v absorption lines, which originate in the atmosphere of the 
primary component, define a rather eccentric binary orbit (e=0.64±0.03). The radial 
velocity variations of the N v absorptions have a lower amplitude than those of the 
He II emission. Such a behaviour of the emission line radial velocities could be due to 
distortions produced by a superimposed absorption component from the companion. 
High resolution cchclle spectra observed during the quadrature phases of the binary 
show H and He n absorptions of both components with a radial velocity difference of 
about 541 km s _1 . From this difference, we infer quite high values of the minimum 
masses, of about 87Mq and 53M Q for the primary and secondary components, respec- 
tively, if the radial velocity variations of the He n emission represent the true orbit of 
the primary. No He I absorption lines are observed in our spectra. Thus, the secondary 
component in the Wack 2134 binary system appears to be an early O type star. From 
the presence of H, He n and N v absorptions, and N iv and C iv emissions, in the 
spectrum of the primary component, it most clearly resembles those of Of/WNLha 
type stars. 
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1 INTRODUCTION 
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tories 
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Inter-American Observatory, Chile, and CASLEO, Argentina 
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iting Astronomer, Las Campanas Observatory and Cerro Tololo 
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f Postdoctoral Fellow CONICET, Argentina 
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** Visiting Astronomer, CASLEO, Argentina 



The X-ray source IE 1024.0-5732, detected with EINSTEIN, 
was identified by Caraveo et al. (1989) with the emission 
line star 2134 in the Wackerling (1970) catalogue, and sug- 
gested to be a binary system composed of an O star with 
a compact companion. Further X-ray data of this source 
obtained with ROSAT were analyzed by Mereghetti et al. 
(1994), who described the optical spectrum to be of Wolf- 
Rayet type, and proposed that the X-rays could arise in the 
colliding winds of a WN+OB binary system. Thus, the star 
was added to the Catalog of Galactic Wolf-Rayet stars (Van 
der Hucht 2001) as WR 21a. Recent interferometric radio 
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observations detected a weak non-thermal source at the po- 
sition of Wack 2134, which was also interpreted as due to a 
colliding wind region in a WN+OB binary (Benaglia et al. 
2005). 

However, indications of orbital binary motion had not 
been found thus far. Here we present a radial velocity study 
of Wack 2134, showing it to be a binary system with an 
orbital period of 31.673±0.002 days and a rather high ec- 
centricity, e=0.64±0.03. 



2 OBSERVATIONS 

We have obtained optical spectroscopy of Wack 2134 be- 
tween 1994 and 2007 at various observatories. A descrip- 
tion of the instrumental configuration of each observation is 
shown in the Table Q] Resolution was determined by mea- 
suring the FWHM of comparison-arc emission lines, and the 
velocity resolution was calculated at 4686 A. 

The spectra obtained at CASLEO, LCO, ESO, and 
CTIO were processed with IRAFjj routines, and those ob- 
tained at SAAO with FIGARO supported by Starlink. 

We have determined the radial velocities (RVs) of 
Wack 2134 by fitting Gaussian profiles to the observed lines 
using the IRAF routine NGAUSSFIT (in the STSDAS pack- 
age). This routine provides an estimation of errors for each 
fitted parameter. For example, in the fitting of the Hell 
4686A emission lines, we obtained errors in line position of 
about 11 km s _1 and 6 km s _1 , for the CASLEO and echelle 
spectra, respectively. 

The spectra were first normalized to the continuum and, 
to standardise the RV measurements in the broad lines, we 
used the position of the core of the line; this procedure has 
the advantage of being less dependent on the errors in the 
definition of the continuum. 

We do not have any target in common among the differ- 
ent instruments, which could be used as a comparison star to 
investigate possible wavelength zero-point differences. The 
different spectral coverages prevented a suitable characteri- 
zation of any such differences by means of interstellar lines. 
However, as will be shown in Section 4, we did not detect 
any significant shifts within our measurement errors. 



3 THE SPECTRUM OF WACK 2134 

Figure [T] shows a low-resolution optical spectrum of 
Wack 2134, obtained at SAAO. Note that C IV emission 
is present in the spectrum, as well as N IV and He II. N ill 
is faint, and N v is observed in absorption. Faint absorption 
lines of H and He n are also present. The observed opti- 
cal spectrum of Wack 2134 resembles that of HD 93162 (= 
WR 25), a WN6ha type star and the second brightest X-ray 
source in the Carina Nebula (NGC3372), which has also 
been recently unveiled as a massive binary system (Gamen 
et al. 2006). The spectrum of Wack 2134 presents intrinsic 
H and He II absorptions (see below) , as well as the C IV and 



2 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
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N IV emissions with comparable intensities. An effort to clas- 
sify this star has been made by Reig (1999). He found spec- 
tral evidence for a WN+Of classification, i.e. a very broad, 
high-intensity He II A4686 emission line but nitrogen (N v 
A4604 and A4620) and hydrogen absorption lines. However, 
the FWHM of the He II A4686 emission line is narrower than 
30 A and, in the following section, we demonstrate that the 
observed absorption lines move in phase with it, so they be- 
long to the same star. Thus, we prefer to classify Wack 2134 
as 03 f*/WN6ha, i.e. as a massive star in an intermediate 
evolutionary stage between an O-type and a Wolf-Rayet star 
similar to those found in the R136 cluster (at the center of 
the 30 Doradus nebula in the Large Magellanic Cloud), cf. 
Melnick (1985), Walborn & Blades (1997), Massey & Hunter 
(1998). Many strong interstellar features are observed in the 
spectrum, indicating that it is a heavily reddened star. 

Spectral lines of the secondary component are observed 
in our high resolution echelle spectra obtained during both 
quadratures of the binary system (see below). These lines 
correspond to absorptions of H and He II. He I lines are not 
visible in our spectra. Therefore, the secondary component 
certainly is an early O-type star, probably as early as 04. 
The spectral type of the secondary may not be much earlier 
than this, since we observe single N v absorptions when H 
and He II appear clearly double. The radial velocity of the 
N v absorptions corresponds to the WN component (as is 
shown below). 



4 THE RADIAL VELOCITY ORBIT 

In all of our spectra we have determined radial velocities 
of the spectral features by fitting gaussians to the line pro- 
files. Only the strongest emission line, He II A4686, could 
be measured in all spectrograms. Weaker features were also 
measured in those spectra with higher S/N. Of these, N v 
absorptions as well as N iv and C iv emissions show radial 
velocities which follow the orbital motion of He II A4686, 
and thus originate in the atmosphere of the primary WN 
component. The radial velocities of the absorption lines of 
H and He II show much scatter, and are separated into two 
components only in our high resolution echelle spectra. The 
radial-velocity measurements of He II A4686 emission, of the 
mean of N v A4604-20 absorptions, and of the N iv A4058 
emission are presented in Table [5] 

We introduced the He II emission line radial velocities 
into a Lafler & Kinmann (1965) period search routine. The 
radial velocities from our spectra do not show large varia- 
tions from one night to the next, but considerable variations 
are present between data obtained during different observing 
runs, thus suggesting a binary period longer than 10 days. 
The best period found was 31.67 days, with some aliases, i.e. 
40.68, 45.64, 63.34 days. An inspection of the distribution of 
the radial velocities phased with each of those periods read- 
ily indicated that the most suitable is 31.67 days and that 
the orbit of Wack 2134 is very eccentric. 

The period of 31.67 days was then introduced as an ini- 
tial value into a routine for defining the orbital elements of 
the binary. To this end we used an improved version of the 
program originally written by Bertiau & Grobben (1969). 
Taking into account the different instrumental configura- 
tions involved in our dataset but also the spectral S/N, we 
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Table 1. Details of the different observing runs. 



Telescope 


Observatory 


Inst. conf. 


Spectral range 

A 


Dispersion 
A pix" 1 


Resolution 
A 


velocity resolution 
km s — 1 


S/N 


n 


1.9-m 


SAAO 


ITS+SITe CCD 


4200-5000 


0.5 


1.5 


96 


16 


6 


1.9-m 


SAAO 


ITS 


3800-7800 


2.3 


6.0 


386 


90 


1 


2.15-m 


CASLEO 


REOSC+Tekl024 


3800-5500 


1.6 


1.2 


269 


80 


24 


2.5-m 


LCO 


Echelle+Tek5 a 


3650-10150 


0.1 


0.14 


8 


10 


8 


4-m 


CTIO 


R-C+Loral3K 


3650-6700 


0.95 


3.8 


213 


100 


13 


3.58-m 


La Silla 


EMMI 


4150-7700 


0.2 


0.6 


38 


80 


1'' 



SAAO: South African Astronomical Observatory; CASLEO: Complejo Astronomico El Leoncito 1 ; LCO: Las Campanas 
Observatory, Chile; CTIO: Cerro Tololo Inter-American Observatory, Chile, 
a: The CCD was binned 2 X 2 to increase the S/N. 
b: Spectrum retrieved from the ESO database. 
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Figure 1. Continuum rectified spectrum of Wack 2134 obtained at SAAO in February 1998. Absorption and emission lines are identified 
below and above the spectrum, respectively. The stellar emission lines identified are the blend of He II with Ha, H/3, H7, and N III+H5: 
N IV A4058, A7103-22; N III A4511-15, A4634, A4640-42; C IV A4658, A5801-12; He II A4686; and the absorptions are N V A4604, A4620; 
He II A5411; and H/3, H7, and H<5. Strong unidentified absorptions are mostly of interstellar origin. 



decided to weight the spectra such that echelle data were 
assigned a value of 1 and the lowest resolution spectrum 
0.1. Thus, SAAO data were weighted with 0.1, CASLEO 
with 0.4, and CTIO with 0.8. We inspected the individual 
O — C values derived by the program looking for systematic 
RV shifts, but in all cases the mean of the O — C (of each 
dataset) remained below 30 km s _1 , which we considered 
as a conservative internal error. Thus we did not apply any 
zero-point corrections to the data. 

We calculated the orbital elements of Wack 2134 inde- 
pendently with the radial velocities of the He n and N IV 
emissions and the N v absorptions. We adopted the period 
determined with the most numerous set of radial velocities 
from the He 11 emission also for the orbit of the N v ab- 
sorptions and N IV A4058 emission. We obtained eccentric 
orbital solutions, e = 0.64±0.03 for the three datasets. Simi- 
lar times of periastron and of maximum radial velocity found 



in each dataset indicate that these lines move together, thus 
belonging to the same component. A similar conclusion is 
reached for the C IV emission lines, confirming that C and 
N lines are formed in the same envelope. The orbital ele- 
ments for the different datasets of radial velocities are listed 
in Table [3] Figure [2] illustrates the radial velocity orbits of 
Wack 2134 as defined by the He 11 A4686 and N IV A4058 
emission lines, and the N v absorptions. In this Figure we 
labeled each instrumental-configuration dataset with a dif- 
ferent symbol in order to show that there are not systematic 
shifts among their RVs. 

We note that the radial velocity orbit defined by the 
He 11 emission has a higher amplitude when compared with 
the radial velocity orbit defined by the N v absorptions. 
The higher amplitude of the radial velocity variations of the 
He 11 emission in principle could arise from distortions of the 
emission line by a superimposed absorption moving in anti- 
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Table 3. Orbital solutions corresponding to the radial velocities of the He II A4686 emission line, the mean of the N V A4604-20 absorption 
lines, the N IV A4058 emission line, and the C IV A5801 emission line in the spectrum of Wack 2134. Symbols have the canonical meanings. 
The last three correspond to, respectively, the mass function, a standard deviation of the fit (parameter defined by Bertiau & Grobben 
1969), and the number of data points involved. 





He II A4686 em 


N V A4604-20 abs 


N iv A4058 em 


C iv A5801 cm 


P[d] 


31.673 ± 0.002 


31.673 (fixed) 


31.673 (fixed) 


31.673 (fixed) 


V [kms" 1 ] 


157 ± 3 


-51 ± 7 


-17 ± 9 


20 ± 5 


K [kms" 1 ] 


172 ± 3 


159 ± 6 


163 ± 8 


138 ± 4 


e 


0.64 ± 0.02 


0.64 ± 0.03 


0.64 ± 0.03 


0.64 ± 0.02 


u> [degrees] 


276 ± 3 


281 ± 6 


287 ± 7 


303 ± 4 


Tperiast [d] * 


190.68 ± 0.08 


191.0 ± 0.2 


191.1 ± 0.2 


191.5 ± 0.1 


TkVmax [d]* 


192.39 ± 0.08 


192.5 ± 0.2 


192.4 ± 0.2 


192.5 ± 0.1 


a sini [Rq] 


82.2 ± 0.5 


76 ± 5 


78 ± 7 


66.2 ± 3 


F(M) [M ] 


7.6 ± 0.1 


6 ± 1 


6.4 ± 2 


4 ± 0.6 


a [kms -1 ] 


16.4 


26.6 


34 


12.5 


n 


52 


36 


10 


23 



*: Heliocentric Julian Date 2,454,000+ 



phase, i.e. a He n absorption originating in the atmosphere 
of the secondary component. However, our spectra do not 
show any clear evidence for an absorption line originating 
in the secondary component. On the other hand, if the sec- 
ondary component were of earlier spectral type, i.e. an 03 
star which also shows N v absorptions in its spectrum, then 
blending of both components could explain the lower am- 
plitude of the radial velocity variations of these absorption 
lines. Spectra of higher S/N and resolution are needed to 
verify these hypotheses. 

Our high-resolution echelle spectra obtained during or- 
bital phases close to the quadratures (<f) ~ 0.93 and 0.05) 
show double absorption lines of H and He II, most clearly 
seen in He n A5411, as illustrated in Figure [3] Measuring 
the component of the He II A5411 absorption line belonging 
to the assumed O companion in the echelle spectra taken 
during quadratures (four spectra obtained during maximum 
RV but only one presenting minimum RV), we obtained a 
difference in radial velocities of about 541 km s _1 between 
quadratures. When we did the same with the absorption 
lines from the WN component, we obtained a radial- velocity 
difference of about 352 km s _1 , which means an orbital semi- 
amplitude of 174 km s _1 , in good agreement with the N iv 
and He n emission lines. Assuming that the radial veloc- 
ity orbit denned by the He II A4686 emission represents the 
orbital motion of the primary (WN) component, and the ra- 
dial velocities of the He n A5411 absorption line (corrected 
by the difference between the systemic velocities of both 
components) show the secondary orbital motion, we per- 
formed a fit of the orbital solution. The new SB2 orbital 
solution implies very high minimum masses for the binary 
components, namely 87±6 Mq for the primary WN type 
component, and 53±4 Mq for the 04: type secondary com- 
ponent. This solution is depicted in Figure [4] If the radial 
velocity orbit of the N v absorptions represents the true 
orbital motion of the primary, then the minimum masses of 
the primary and secondary components would be 83±22 Mq 
and 47±14 Mq, respectively. Similar values were obtained 
when we used both components of the He II A5411 absorp- 
tion line as representing the orbital motion of each star in 
the system. 
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Figure 3. He II A5411 absorption lines of both primary (WN) 
and secondary (O) components in the Wack 2134 binary system 
observed in our high resolution echelle spectra during the quadra- 
ture phases of the orbital motion. We also show (in the middle) 
a spectrum taken at nearly conjunction when the lines from both 
components are blended. Note that the radial velocities are more 
extreme in the O-type star, thus indicating a lower mass. Diffuse 
interstellar bands (DIB), which are blended with the stellar lines 
in lower resolution spectra, are also indicated. 



With the high minimum masses found from our ra- 
dial velocity orbit, we would expect to observe eclipses. 
Wack 2134 has been monitored in V magnitude by the All 
Sky Automated Survey (AS AS) (cf. Pojmafiski 2001). We 
have examined the public ASAS data of this star. No ob- 
vious variations were found in the V magnitudes when we 
folded them at the spectroscopic binary period of 31.673 
days. However, as Wack 2134 is near the faint magnitude 
limit of the ASAS survey, the V magnitudes show rather 
high noise. More accurate photometry is needed to rule out, 
or confirm, the eclipsing nature of this binary. 
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Figure 2. Radial velocity variations of the He II A4686 emission, N V A4604-20 absorption, Niv A4058 emission, and Civ A5801 emission 
lines in the Wack 2134 binary system, phased with the period of 31.673 days. We distinguished the RVs measured in the different datasets: 
SAAO spectra (x), EMMI spectrum (•), CASLEO (o), CTIO (■), and LCO (A). Note that there are no systematic shifts among the 
data obtained with different instrumental configurations. The continuous curves represent the orbital solutions from Table [3] 



5 CONCLUSIONS 



From the radial- velocity variations of the He II A4686 emis- 
sion line in the spectrum of the WN-type star Wack 2134, 
we found a most probable period P=31.673±0.002 days, 



thus revealing Wack 2134 as an eccentric binary system 
(e=0.64±0.03). 

Our higher-resolution spectra taken during quadratures 
show H and He II resolved into two components. We could 
obtain an orbital solution for the secondary component mea- 
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Table 2. Radial velocities of some lines in the spectrum of Wack 
2134. 



HJD a 


He II 


N v 


N iv 


Civ 




A4686 em 


A4604-20 abs 


A4058 em 


A5801 em 




[km s _1 ] 


[km s- 1 ] 


[km s- 1 ] 


[km s- 1 } 



0179.427 


272 








0621.234 


248 








0623.269 


213 


35 






0850.557 


129 






18 


0879.420 


184 


-22 






1189.580 


323 








1484.854 


139 


-90 






1653.528 


3 


-189 


-45 




1654.516 


4 


-171 


-68 




1655.519 


27 


-211 


-78 




2009.574 


282 




130 




2011.589 


285 




47 




2013.561 


210 


20 


75 




2298.747 


251 








2299.688 


191 


-61 


158 




2353.527 


161 


17 




54 


2384.512 


53 




-79 




3076.798 


33 


-228 


-73 




3077.790 


20 




-88 




3145.540 


202 


-128 


11 




3146.483 


282 


15 


186 




3150.519 


246 


29 







3151.476 


246 




40 




3154.477 


163 


-39 


85 




3155.463 


185 


-2-1 


-4 




3156.479 


187 




52 




3169.453 


64 








3170.449 











3171.453 


30 


-189 


-215 




3172.459 


32 




-263 




3481.536 


115 


-66 


-62 


-14 


3482.538 


119 






-57 


3489.501 


28 


-132 


-117 


-68 


3490.540 


IS 


-156 


-143 


-85 


3491.540 


20 


-160 


-143 


-71 


3741.788 


3 


-179 


-153 




3747.790 


325 




82 




3772.701 


21 


-125 


-112 


-88 


3875.510 


311 


188 


255 


196 


4188.714 


-11 


-209 


-200 


-18 


4188.724 


-8 


-216 


-160 


-32 


4189.508 


11 


-209 


-52 


-28 


4189.517 


49 


-167 


-194 


-15 


4189.769 


-1 


-187 


-156 


-50 


4189.778 


11 


-196 


-184 


-53 


4190.525 


124 


-129 


-47 


21 


4190.801 


258 


-16 


12 


47 


4191.510 


297 


11 


76 


116 


4192.489 


3 15 


126 


179 


202 


4192.765 


3 15 


124 


130 


211 


4193.479 


363 


122 


165 


193 


4193.794 


332 


121 


99 


175 


4200.583 


226 


8 


61 


60 



a: Heliocentric Julian Day: 2,450,000+ 



























_ 


~ 400 




* 
















\ , 




- 




* * 


'- 




= 












\ *£~* 


= 




c 300 






-"■ = 
















•\= 


- 


X 


_ 




M^ 


"~5 ^ 














s^i% 


- 


>, 200 








^%w 


\S 












W^t 


ft — 









































^ 


;*. 












■f, 100 




















\ 






> 


















k~ s 






- 


| 


J" 














•V 




1 1 




- 


(0 


























K-100 




















i 




- 



0.5 1 

Orbital Phase 

Figure 4. Radial velocity variations of the He II A4686 emission 
and He II A5411 absorption line, in the Wack 2134 binary system, 
phased with the period of 31.673 days. Filled symbols represent 
the radial velocities of the emission (WN component), and open 
ones depict the RVs of the absorption (O component). Curves 
show the orbital motion of each component. 



suring the He II A5411 absorption line in those spectra, but 
this orbit has only one measurement at one of the quadra- 
tures (see Fig. [4}. Using the He II A4686 emission line and 
the He II A5411 absorption line to follow the motion of the 
WN and O component, respectively, we obtained very high 
minimum masses Mwn~87M0 and Mo~53M0. In spite of 
the high minimum masses obtained, no light variations in- 
dicative of eclipses are found in the photometry performed 
by AS AS. These minimum masses have to be considered as 
preliminary, as more high resolution and S/N spectra are 
needed to confirm these values. Also, more precise photom- 
etry could shed light on this new massive binary system. 
Although we prefer to be cautious, the WN component in 
the Wack 2134 system could be amongst the most massive 
stars ever measured from binary motion in the Galaxy, as 
both WN6ha components in the binary system WR 20a 
(83+82M ; Rauw et al., 2004; Bonanos et al., 2004), lo- 
cated about 16' from Wack 2134, and those in the new 
Wolf-Rayet binary system NGC3603-A1, which has masses 
of 116+89M Q (Schnurr et al., 2008; Moffat et al., 2004). 
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